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Abstract 

We review recent results obtained for the dynamics of incipient chaos. These results 
suggest a common picture underlying the three universal routes to chaos displayed 
by the prototypical logistic and circle maps. Namely, the period doubling, inter- 
mittency, and quasiperiodicity routes. In these situations the dynamical behavior 
is exactly describable through infinite families of Tsallis' (^-exponential functions. 
Furthermore, the addition of a noise perturbation to the dynamics at the onset 
of chaos of the logistic map allows to establish parallels with the behavior of su- 
percooled liquids close to glass formation. Specifically, the occurrence of two-step 
relaxation, aging with its characteristic scaling property, and subdiffusion and arrest 
is corroborated for such a system. 

Key words: Nonlinear dynamics, Renormalization group. Weak chaos, Glass 
formation 
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1 Introduction 

In recent years there has been a renewed interest in the analysis of the dynam- 
ical behavior at critical attractors in low dimensional maps. In such situations, 
criticality is triggered by the vanishing of the Lyapunov coefficient (coefficients 
in dimension higher than one) so that ordinary chaos paradigms are not ap- 
plicable. Part of this interest is because Tsallis' functional forms (the so-called 
g-exponential and g-Iogarithm) were found to mimic the role that the ordi- 
nary exponential and logarithm have in the usual theory of chaos (although, 
as we will see below, in a more complex and sophisticated sense). In fact. 
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exact renormalization group (RG) solutions for the sensitivity to initial condi- 
tions at these critical attractors are found to assume the structure of (infinite 
families of) g-exponential functions. This paper attempts to rationalize in a 
single review some of the most relevant new results that stemmed from the 
application of exact RG approaches to critical attractors of low dimensional 
maps [1,2,3,4,5,6,7,8]. This research work moved from the understanding that 
RG techniques developed in the '80s for describing the static properties of the 
maps could indeed be extended in order to solve in all details the dynamical 
behavior at critical attractors. Then, contact with TsaUis' and Mori's for- 
malisms were established a posteriori. The final result is the identification of a 
universal, intricate, dynamical mechanism common to the three paradigmatic 
routes to chaos of period doubling, intermittency and quasiperiodicity. 

Another important result conveyed by the RG approach at the chaos threshold 
of the logistic map is the determination of the exponent characterizing the 
decay of correlations in presence of an external noise source [7,8]. By realizing 
that such an iteration equation is a discrete version of Langevin's equation and 
that the underlying scenario is that of an ergodic to non-ergodic transition, one 
can in fact produce a "translation dictionary" between the dynamical behavior 
of the logistic map with additive noise and the phenomenology observed for a 
glass former. Translated words include Adam-Gibbs formula, time translation 
invariance and a-relaxation, aging, subdiffusion and arrest. 

There is a further interesting mechanism through which aging and weak er- 
godicity breaking can be obtained in low dimensional maps. This mechanism 

is related to the dynamical behavior close to the origin of a Manneville-like 
map (see below) and can be statistically described in terms of the continuous 
time random walk formalism. Such a discussion is however out of the scope of 
the present paper. The interested Reader is referred to [9] for an appropriate 
account. 

The paper is divided in two main sections. Section 2 is devoted to the dynam- 
ical picture underlying the above mentioned routes to chaos. In Section 3 we 
describe the analogies between the behavior of the logistic map with additive 
noise and that of a glass former. Finally, some remarks are added. 



2 Routes to chaos: a common dynamical picture 

In the following, we will be mainly considering unimodal maps, i.e., discrete- 
time iterated maps governed by the equation 

xt+i^ ff,,Q{xt) = i- n\xt\^, f = 0, 1,..., xe[-l, 1], (1) 
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where // e [0, 2] is the control parameter and C > 1 the non hnearity order of 
the maximum at x = (Fig. 1). C, = 2 corresponds to the celebrated logistic 
map [10,11]. We define the sensitivity to initial conditions as 

i{xQ,t)= lim 1^^, (2) 



where we denote explicitly the dependence on the initial condition xq in order 
to address cases in which the Lyapunov coefficient is zero. Typically (for all 
II except a set of zero Lcbcsguc measure), ^ is independent from xq and ex- 
ponential for t » 1, ^{t) ~ exp (At), where A e M is the Lyapunov coefficient 
[10,11]. 

For smaller than a critical value fiociC), ^ is negative except for an infinite 
but numerable number of pitchfork bifurcations points in which A = (see Fig. 
2). Within this region, the attractor is periodic of order 2^ (r = 0, 1, . . .) with 
r that increases by one unit at each pitchfork bifurcation and that becomes 
infinite at ^oo- Inside each //-interval in which the period is 2'', there is a 
special orbit that passes through a; = and that has A ^ — oo (super-stable 
cycle [10]). Important quantities for the dynamics at the onset of chaos are 
the diameters of the bifurcation forks dr,s {s — 0, 1, . . . , 2'' — 1), which are 
defined as the distances of two nearest neighbors attractor points in a super- 
stable cycle. As we will see, in the limit r ^ oo these diameters contain all 
the information about the multifractal attractor at the onset of chaos. The 
chaos threshold is the accumulation point of both pitchfork bifurcations 
and super-stable cycles. 

For n > Hoc, A is positive except in an numerable number of intervals which are 
again characterized by a periodic attractor. The left extreme of such intervals 
is determined by a tangent bifurcation where A vanishes. Inside each interval, 
a new pitchfork bifurcation cascades produces a (higher order) chaos threshold 
before reaching the right extreme of the interval, so that the final result is an 
intricate and fascinating self-similar structure. 

The study of unimodal maps marked an era in chaos theory. For example, the 
pitchfork cascade has been identified to occur in chemical reactions, optically 
bistable systems, electrical RCL oscillators, sound waves in water, Bcnard con- 
vection fiuid experiments, . . . [10]. Also, the tangent bifurcation with the addi- 
tion of a reinjection mechanism is the paradigm for the phenomenon of (type 
I) intermittency and l//-noise observed, e.g., in nonlinear RCL oscillators and 
in Benard convection [10]. The elegance and richness of this self-similar struc- 
ture stimulated a number of theoretical results starting from the late 70's. 
Perhaps the milestone was the RG approach developed by Feigenbaum [12] 
and, independently, by CouUet and Tresser [13] for the pitchfork bifurcations 
cascade. This work was extended by the group of Pohti [14] and that of Mori 
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Fig. 1. Unimodal map. 

[15] that focused on the properties of the fluctuation spectrum of generahzed 
algebraic Lyapunov coefficients. An exact solution for the tangent bifurcation 
dynamics was found by liu and Rudnick [17]. As we show below, such solu- 
tion can be extended also to the pitchfork bifurcations. Another important 
contribution by Mori and colleagues was the introduction of the idea of dy- 
namical transitions [16], a formalism closely related to the probabilistic large 
deviation theory [18], that is analogous to the description of first-order ther- 
mal transitions. More recently, the group of Tsallis [19] pointed out that the 
g-exponentials derived within the so-called nonextensive statistical mechanics 
formalism appropriately describe the envelope of the fluctuating sensitivity 
to initial conditions in which the ordinary Lyapunov coefficient vanishes. The 
work by Roblcdo and colleagues [1,2,3,4,5,6] that we are reviewing in this sec- 
tion put together all previous approaches in a single perspective, by deriving 
exact results supported by the KG functional composition. 



2.1 Pitchfork and tangent bifurcations 

Let us start by understanding what happens at pitchfork and tangent bifur- 
cations. Reading such a transition along the /i-axis of Fig. 2a, the transition is 
said of order n if in the left (pitchfork) or in the right (tangent) neighborhood 
of the transition point the attractor is periodic of order n. If we consider the 
n-composed map f^^^ and shift the coordinate x to one of the of the n bifur- 
cation points ^ , we have that in the neighborhood of the chosen point f^^^ has 



With an abuse of notations we still indicate the shifted coordinate X ctS X — X. 
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Fig. 2. Attractor (a) and Lyapunov structure (b) of the logistic map (C = 2). 
the universal expansion 

/S(^) = ^ + (-1)^^ sgn{x^)\x\^ + o{\xn, (3) 

where w > is the leading expansion coefficient and sgn is the sign function 
(here and below repeated use of the sign function is made in order to unify the 
treatment of pitchfork and tangent bifurcations). For a pitchfork bifurcation 
z = 3, whereas for a tangent bifurcation z = 2 (see Fig. 3). The expansion is 
independent of the non-linearity parameter (. 

With z — 2, hy neglecting the o(|a;|^) term, Eq. (3) recovers the Manneville 
map [10]. In such at the left side of the bifurcation point trajectories 

are converging to the bifurcation point, whereas at the right side they are 
expelled from the bifurcation point (Fig. 3b). As we will see below, in the 
converging side the dynamics is regular, being characterized by a weak (power 
law in place of exponential) insensitivity to initial conditions. In contrast, the 
expelling side is very strongly divergent, with a sensitivity to initial condi- 
tions stronger than exponential. Thus, the addition of an external reinjection 
mechanism connecting the right hand side with the left hand side allows an in- 
termittent dynamics in which chaotic and regular patterns alternate. Recently, 
a remarkable connection between this mechanism and the dynamics of a clus- 
ter of spins at a thermal critical point has been established [20]. This opened 
the possibility for the description of such temporary critical fluctuations in 
terms of Tsallis' functional forms [21]. 

In correspondence with the tangent bifurcation (z — 2), Hu and Rudnick [17] 
discovered that the RG iteration equation 
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Fig. 3. Pitchfork (a) and tangent (b) bifurcations. 

possesses a solution, f*{x) = x[l — (z — l)u sgn(x)\x\^~^]~^^''^~^\ that coincides 
with the expansion (3) up to order o(|x|^). Remarkably, this is one of the very 
few RG solution exactly known [11]. The original solution [17] can be extended 
also to the pitchfork bifurcations by simply including a change of sign [2]: 

f*{x) = x[l -{z-l)u {-ly sgn(x") sgn{x)\x\'-Y'/^'~'\ (5) 



Close to the bifurcation point, by using the property: 

(^) = ^TT^rTyr K^^^"'^^), m = 1, 2 . . . , (6) 



m 



we can recast the m-times iterate of the n-composed map /^"^ as 

(m) 



xo[l -{z- l)a{-iysgn{x^,)sgn{xo)\xorH]-'/^^-'\ (7) 



where t = mn and a = u/n. Notice that while the n-composed map describes 
a single orbit of this form in the shifted coordinate xq <^ 1, for the original 
map ffj_^c there are n specific time subsequences, each corresponding to one of 
these orbits originated in one of the n bifurcation points. Eq. (7) satisfies the 
identity 

dxt _ / xt_ 
dxo \xo 
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Fig. 4. g-exponential function for positive and negative arguments and different 
values of the parameter g G R. 

Due to this property the sensitivity to initial conditions ^(xq, t) = Huiaxo^o 
assumes the form 

e(xo,t) = [l + (l-g)A,(xo)t]^ =expg(A,t), g e R. (9) 



The parameter q can be read from z. For all pitchfork (tangent) bifurcations 
of all order n we have q = 5/3 (g = 3/2). The functional form exp^ is a one- 
parameter deformation of the ordinary exponential function called Tsallis' 
g-exponential [19]. In the limit g — 1 the ordinary exponential is recovered. 
The inverse of the g-exponential, named g-logarithm, is given by lng{x) = 
{x^~'^ — l)/{l — q). The coefficient Xg replaces the ordinary Lyapunov coefficient 
and depends on the shifted coordinate of the initial point: 

Xg{xo) = za {-ly sgn(xg) sgn(xo)|xo|^"^ (10) 



As for the ordinary Lyapunov coefficient, the sign of Xq discriminates if the 
dynamics is regular {Xq < 0, insensitivity to initial conditions, meaning that 
two different orbits tend to converge in time) or chaotic (A^ > 0, sensitivity to 
initial conditions, i.e., two different orbits tend to separate in time). In addi- 
tion, the value of g determines whether the dependence with time is stronger 
or weaker than exponential. Fig. 4 reports the behavior of the g-exponential 
for different g's and sign of A^. 

In the next subsection we will see that a similar scenario characterizes the 
onset of chaos. Namely, the ordinary exponential form for the sensitivity to 
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Fig. 5. Time subsequences generated at the onset of chaos from x-j-^q = 0. Numbers 
label the iteration time r. Each subsequence lies along a different diagonal line. 

initial conditions, which is typically independent of the initial condition Xq, 
is replaced by an infinite family of Tsallis g-exponential functions, each as- 
sociated to a specific time subsequence, and where the generalized Lyapunov 
coefficients Xq become dependent on xq. Such intricate dynamical picture is 
then common to all pitchfork and tangent bifurcations, and to the periodic 
and quasi-periodic onset of chaos. 



2.2 Onset of chaos 



For the sake of simphcity we address first the logistic map C = 2. We also 
change slightly our notations indicating for a moment with r = 0, 1, . . . the 
iteration time. Eq. (1) reads then 

Xr+l = 1 - fJ^xl- (11) 



The chaos threshold is located at /ioo(2) — 1.40115 Let us start by choosing 

the initial condition xq — 0. The Feigenbaum attractor can then be character- 
ized [3,4] by a series of monotonic subsequences 

Tfe = (2A; + l)2"-^ A; = 0,1,..., n>k. (12) 



The log- log plot of Fig. 5 shows that for each fixed k, n > k defines a subse- 
quence such that the absolute values of the iterates \xri, \ lie along a diagonal 
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line. 



The Fcigcnbaum-CouUet-Tresser RG transformation defines the universality 
class of the onset of chaos for the logistic map. Such transformation is given 
by [10,11] 

g{x) = a''g^^"\x/a''), a = 2.50290 ... , (13) 



where a is one of the Feigcnbaum's universal constant and g is the universal 
Feigenbaum's function (not known in closed form). This RG transformation 
allows to calculate explicitly the iterates for each subsequence r^. In fact, one 
obtains [4] the expansion 

. , o(2fc+l)CQ^ (2fc+i)//CQ\ 

-^.-IS'^'WI-^J + ^^J^. (14) 



which is valid for initial conditions |a;o| < a We now shift to zero the 

initial time of each subsequence by redefining the time as 

tk = Tk - 2k - 1 = {2k + 1)2''-'' - 2k - 1, n>k. (15) 



Implementing this time shift in Eq. (14) and using the identity a^""^) = 
[1 + tk/ {2k + i)]in«/in2^ ^Yie sensitivity to initial conditions for each time sub- 
sequence turns out to be [4] 



exp. 
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In a 
(2A;+ l)ln2 



(16) 



(see Fig. 6). The fact that the generalized Lyapunov coefficient Af varies with 
the subsequence k can be seen as a dependence of this coefficient on the initial 
condition. In fact, for each subsequence, the initial position at zero time tk, 
Xr=2k+i, changes with k. 

Along these lines, it is also possible to show [4] that within each subsequence, 
Tsallis' g-entropy grows linearly with the time tk and the slope of the linear 
growth is equal to ^ Af . This, parallels the Pesin identity that relates the (sum 
of the) positive Lyapunov coefRcient(s) to the Kolmogorov- Sinai entropy rate 
[11], and what happens for the ordinary entropy growth {q — 1) when strong 
chaos is present. 

^ The reader may be interested by the fact that in a different context, namely 
for a conservative bidiniensional map with vanishing Lyapunov coefficients, similar 
numerical evidence has been found, although for the somewhat trivial (linear) case 
q = [22]. 
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Fig. 6. Sensitivity to initial condition at the onset of chaos of the logistic map. 
For different subsequences, the g-logarithm/linear plot with g = 1 — exhibits 
straight lines whose slope corresponds to the generalized Lyapunov coefficient \q. 



Although Eq. (16) gives the sensitivity to initial conditions for all the subse- 
quences A;, it is not yet completely general since such subsequences are those 
generated starting at a;^=o = 0. At this point, we have to say what happens if 
we choose another point of the attractor as the starting point. The answer is 
given by considering the Feigenbaum's a function. This function is defined by 
the scaling properties of the diameters of the bifurcation forks dr^s'- 



ar{s) 



d 



(17) 



In the limit r oo. a,, has an infinite number of jumps, each associated 
to a different attractor point which is here identified by the variable s. The 
sensitivity to initial conditions for the subsequence k generated by the initial 
point Xt=o defined by s can be deduced from a through [5] 



lim 

n—*oo 



dn-k+l,s+t 




CT„_fe(s - 1) 


dn—k,s 




CTn-kis) 



n—k 



n—k 



(18) 



where A is a number that depends on s (and therefore on Xt-=o)- The quantity 
A has been identified as the ratio of the values of a at each of its discontinuities 
[5]. Hence, again the use of the identity A^ — 2^X1]^ lead to the more 
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general result 



ln2 InA 

h^' (2A; + l)ln2- ^ ' 



We see then that the structure valid for Xr=o = is maintained, but changing 
the starting points implies a change of the power law exponent q. 

Finally, all previous results for the onset of chaos can be further generalized 
to general nonlinearity ^ in unimodal maps [5], in which case the families of 
g-cxponential functions characterizing the sensitivity to initial conditions are 
defined by 

^ - ^ " (C-l)lnA' ^'^ - (2A; + l)ln2' 



where A now also depends on 

Summarizing, we have seen that the sensitivity to initial conditions at the 
chaos threshold can be understood in terms of a complex hierarchical struc- 
ture that can be exactly deduced from the Feigenbaum-Coullct-Trcsscr RG 
approach. For a given point of the attractor, this structure corresponds to an 
infinite family of g-exponentials, each with the same q but with a different 
generalized Lyapunov coefficient A^*^^ . Spanning the initial point along the at- 
tractor corresponds to an infinite number of such families, each of the families 
being now characterized by a different value of q which is given in terms of 
the discontinuities of Feigenbaum a function. Since the stronger of such dis- 
continuities happens for s = {xr=o = 0), the value of q corresponding to this 
discontinuity is in some sense the most important one. Such a value, for the 
logistic map, is = 1 — In 2/ Ina = 0.2445 

Another fundamental result [5] for the period doubling onset of chaos (and 
actually also for the quasi periodic route to chaos [6]) is that the TsalUs' q- 
exponential structure for the dynamics can be exactly linked with the formal- 
ism of Mori's transitions [16], which we briefly sketch below. Mori's formalism, 
when applied to the transition to chaos, starts by defining a different gener- 
alization of the Lyapunov coefficient, called generalized finite-time Lyapunov 
coefficient 

1 



(21) 
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The probability for getting the value A at the finite-time t is postulated in the 
form 

p{x,t)^r^^Mo,t). (22) 



This allows the introduction of a dynamic partition function Z{q, t) and of 
the associated free energy 0(q) 



Z{q, t)^ f dXp{X, t)t-^^-'^\ 0(q) ^ - hm i^i^M, (23) 
J t->oo In t 



where q e R is Mori's q-parameter. The "coarse grained" function of gener- 
alized finite-time Lyapunov coefficients is 

A(q) . (24) 



and ip{X) is the Legendre transform of 0(q). In analogy with ordinary ther- 
mal phase transitions, Mori and colleagues detected a single first order phase 
transition for the onset of chaos [16]. 

The important identification [23] 

A = — In A"-'^ = - ln„ (25) 
Inn n 



together with the previous results, leads to the conclusion [5] that the dynamics 
at the onset of chaos is actually characterized by an infinite family of Mori's 
q-phase transitions each occurring at one of the allowed values of the Tsallis' 
q-mdex in the sensitivity: 

q = ?- (26) 
2.3 Quasiperiodic onset of chaos 



The prototype for the quasiperiodic route to chaos is the circle map (see, e.g., 
[24]) 

Ot+i = fnM^t) ^et + n-^ sin(27r^t) (6 mod 1), (27) 

ZTT 



which depends on the real parameters fl (bare winding number) and K (amount 
of nonlinearity) . In order to describe the dynamical behavior it is important to 
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consider the influence of the parameter Q through the dressed winding number 
w ^ hm (28) 



that gives the average increment of 6t per iteration. The orbits for K < 1 
are periodic (locked motion) if w is rational and quasiperiodic (unlocked mo- 
tion) if w is irrational. In fact, w(f2) describes the so-called "devil staircase" 
making a step at each rational value of Q [24]. K — 1 defines the critical 
circle map for which the periodic motion covers the whole Q interval apart 
from a multifractal set of values for which the motion is quasiperiodic. For 
K > 1 regions of periodic motion (Arnold tongues) overlap in chaotic bands. 
The chaos threshold can be studied by fixing K — 1 and by approximating a 
target irrational value of w through a sequence or rational values. Perhaps the 
most interesting case is the one in which the target value is the reciprocal of 
the golden mean Wgm = {V5- l)/2 and the sequence of rationals is given by 
Fibonacci numbers Wn = Fn-i/Fn (with = Fn + F„_i) [24]. The univer- 
sal properties of the onset of chaos are then described by the fixed-point map 
ggm{G) of the following RG transformation 

CX.gm 9gm I Oigmdgm 

{-^ I (29) 

\ \^gmj / 



where agm — —1.288575 ... is a universal constant. The same techniques im- 
plemented for the pitchfork bifurcation onset of chaos lead to the previous 
universal picture, although this transition to chaos is much more involved 
and reach in detail [6]. Namely, the sensitivity to initial condition at the chaos 
threshold is given by an infinite number of Tsallis' g-exponential families, each 
g-exponential being associated to a specific time subsequence ^ : 



itk = exp 



q = l + 



\nWgrn_ 

21nUr 



q 



' k In Wgm ' 



(30) 



where A depends on the initial condition that generates the time subsequence 
k {A — if the initial condition is at = 0). Correspondingly, the dynamics 
is characterized by an infinite family of Mori's q-transitions occurring at q — q 
[6]. 



^ More precisely, the time subsequences are in this case identified by two different 
numbers k = 1,2,... and I = 0,1, . . . , k — 1. Eq. (30) is valid if I = (see [6] for 
details) . 
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2-4 Tsallis' q-exponential sensitivity to initial conditions 

We have seen that if one is interested in extracting significant trends for the dy- 
namics at critical attractors where the ordinary Lyapunov coefficient vanishes, 
some precise time subsequences that depend on the initial condition must be 
selected. Within these subsequences, Tsallis' g-exponentials and g-gcncralized 
Lyapunov coefficients Xg appear to be convenient tools for extracting these 
trends. Table 1 summarizes the different behaviors obtained from various com- 
binations of q and A^. The only case not yet observed is < 1 and Xg < 0, 
which would amount to super-strong (stronger than exponential) insensitivity 
to initial conditions. In private conversations with A. Robledo, we conjectured 
that this could be the case of the super-stable cycles of unimodal maps. Efforts 
to clarify this point are currently in course [25] . 





9<1 


q = l 


q>l 


Ag < 


super-strong 
insensitivity 
super-stable cycles? 


strong 
insensitivity 

IJ'< fJ-oo 


weak 
insensitivity 
pitchfork, tangent Ihs 


Xg>0 


weak 
sensitivity 
M = y"oo 


strong 
sensitivity 
typically for n> Hoo 


super-strong 
sensitivity 
tangent rhs 



Table 1 



Sensitivity to initial conditions and Tsallis' g-exponentials. 

3 Logistic map with additive noise and glassy dynamics 

The theoretical understanding of the dynamical and thcrmodynamical fea- 
tures of supercooled liquids is nowadays a very important research topic in 
condensed matter physics and in statistical mechanics. Theoretical efforts are 
motivated by the search for general common features of non-equilibrium sys- 
tems and applications span from new material science to protein folding. In 
a glass former there are two different dynamical scenarios, separated by the 
glass transition temperature Tg. Above Tg the system equilibrates while below 
Tg it is out-of-equilibrium. Phenomenology of glassy systems includes [26]: 

• Time translation invariant relaxation processes. Above, but close to, the 
glass transition two-time correlations do not depend on the first (waiting) 
time and display power law decays (a-decays). Some glass formers also 
exhibit a further initial power law decay called /3-decay. 

• Loss of time translation invariance (aging) at and below the glass transition. 
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Intriguing connection between kinetics and thermodynamics - Adam-Gibbs 
formula: 



tx — A exp 



( 



B 



) 



(31) 



where t.^,. is a relaxation time (equivalently the viscosity), A and B are 
constants, T is temperature, and Sc configurational entropy 
• Critical slowing down: transition from normal- to sub-diffusion to localiza- 
tion of molecular motion. 

In order to elaborate the counterpart of such behaviors in low- dimensional 
maps, let us start by noticing that in the Langevin theory the effect of collisions 
of the diffusive particle with molecules in the fluid is represented by an additive 
noise term. In the same spirit, nonlinear low-dimensional maps with external 
noise can be used to model systems with many degrees of freedom. The discrete 
form for a Langevin equation is 

Xr+l ^ Xr + h^{Xr) + (JFt- T = 0, 1, . . . , (32) 

where is a Gaussian white noise ((r,-r^/) = S^^') and a measures the noise 
intensity. Taking h^{x) — l — x — jix^, Eq. (32) becomes the iteration equation 
of the Logistic map with additive noise 



For small noise amphtudes the onset of chaos )Uoo(c) still separates an "equihb- 
rium" ergodic phase characterized by a chaotic dynamics with positive Lya- 
punov coefficient from an "out of equilibrium" non-ergodic one where the 
Lyapunov coefficient is negative [10]. 

We will be studying the chaos threshold from /j, > /j,^. In our analogy with 
the glass former this means that we will be looking at the glass formation 
from the side of the liquid. In the absence of noise, the attractor, besides 
being the accumulation point of the bifurcations cascade, can also be viewed 
as the accumulation point of the band splittings occurring above the onset 
of chaos. For — 2, the chaotic attractor is in fact characterized by a single 
band that spans the interval [—1,1] (see Fig. 7a). As /i is reduced towards 
this unique band reduces its width and then splits in two disjoint bands. This 
band splitting continues up to ii^o) where the attractor can be viewed as 2°° 
disjoint bands. The same occurs if cr 7^ 0, the only difference being that the 
band splittings end (and the bifurcations end) at a finite value 2-^*^'^) (Fig. 7b). 
This phenomenon is called bifurcation gap. 

If the noise term is small enough, the RG perturbative approach can be applied 
in order to calculate the position of the iterates around Xq — Q [7] . The result 



Xr+i — 1 — iJ.x1-\- aVr^ xe[— 1,1], /xG[0,2]. 



(33) 
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Fig. 7. Bifurcation gap phenomenon, (b) and (c) are magnification of the region 

enclosed by the box in (a). The infinite number of bands (points) at the noiseless 
chaos threshold in (b) become a finite number of chaotic bands in presence of the 
noise in (c). 



IS 



1-9 



(34) 



where |a;o| <C 1, G\{x) is the first order perturbation eigenf unction, r ~ 0.6332 
is determined from the noisy power spectrum of the map [10] , and we are using 
a notation consistent with the one in Section 2.2. For the leading subsequence 
/c = 0, by introducing the previous time shift t = to = tq — 1, the iterate 
positions can be viewed in terms of g-exponentials: 



xt = exp2_„ {-Xqt) [l+o-rjexp^ (-Kt)] ■ 



(35) 



This implies the presence of a "crossover" or "relaxation" time in the noisy 
dynamical behavior: 

= a'-\ (36) 



For t < tx — 2^^'^'> — 1 the dynamics is restricted to non-overlapping bands 
attr actor, similarly to the noiseless case; for t > the dynamics becomes 
more and more chaotic as a consequence of the band merging (Fig. 8). 

This crossover furnishes the key for identifying the following marks of glassy- 
like behavior for the noisy onset of chaos (see [7,8] for details). 
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Fig. 8. Crossover dynamics observed for the main subsequence A; = for two 
different values of the noise intensity. For t < tx{(j) the iterates closely follow the 
noiseless attractor. For t > tx{(T) the iterates become more chaotic and they jump 
among different bands of the noiseless attractor. 

• Analogue of the Adam-Gibbs formula. The entropy associated to the distri- 
bution of the iterate positions within the 2^ bands has the form 



where as = —a J dV p(r) Inp(r) is the entropy associated to a single band. 
Since 2^ = 1 + and = cr^'-^) we get [7]: 



Notice that in contrast with the Adam-Gibbs exponential law, in this case 
tx has a power law behavior in 1/Sc- 

• Time translation invariance and a-relaxation. Whereas a study of the re- 
laxation processes of iterates starting with initial conditions Xq outside the 
attractor could furnish some analogue to the f3 initial fast relaxation ob- 
served for some glass formers, the a subsequent relaxation is mimicked by 
the bifurcation-gap crossover for t > t^. This becomes apparent if one stud- 
ies the behavior of ensemble-averaged correlations Ce{ti,t2) [8]. There is a 
"thermalization" process occurring for small values of ti. If ti is large enough 
inside the chosen subsequence k, the correlation becomes time translation 
invariant and characterizes the decay of correlations. 

• Aging. At the onset of chaos when a = trajectories are non-ergodic and 
retain memory of the initial data. This property is equivalent to a loss of 




(37) 



l-r 




(38) 
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time translation invariance and to a "built-in" aging. Since ensemble and 
time averages are not equivalent, a time- averaged definition of correlations 
allows the exhibition of such a property [8] . 

Subdijfusion and arrest. The sharp slowing down of dynamics in supercooled 
liquids implies a crossover from normal diffusion to sub-diffusion and finally 
to arrest. This deceleration is caused by the confinement of any molecule 
by a "cage" formed by its neighbors. Such a "cage" can be reproduced by 
a periodic map obtained via repetition of a single cell map in such a way 
that diffusion is due only to the noise term. The periodic map is given by 

Xt+i^F{xt), F{l + x)^l + F{x), ; = 1,0,1,..., (39) 

where 



F{x) 



|1 — llc^"^] + (7^, —1<X< 0, 

|1 - HcX^l +a^, <x <1. 



(40) 



The diffusion process displayed by the map in Eq. (39) reproduces the 
crossover normal/sub-diffusion/ arrest as a tends to zero [8]. 



4 Final remarks 



Several nonlinear dissipative systems, like those mentioned in this review, 
exhibit both regular and chaotic behavior depending on some external param- 
eter. The transition from regularity to chaoticity implies a dramatic alteration 
in the dynamical behavior, the situation changing from an exponentially fast 
insensitivity to initial conditions to an exponentially fast sensitivity to initial 
conditions. There are only few possible mechanisms through which such transi- 
tions occur. These are the well known (universal) routes to chaos. Thus, when 
the external parameters are tuned at the transition point, the ordinary Lya- 
punov coefficients vanish and they do not allow to extract significant trends 
characterizing the critical dynamics. Under these conditions, intricate oscilla- 
tions typically appear. Since these oscillations are reminiscent of the regular 
behavior, the possibility arises to characterize the dynamics exactly if one 
accepts to describe it via a set of specific temporal subsequences. 

Here we have given an account of the extension of the static RG approaches, 
discovered in the 80's, to the dynamics at the critical attractors of period 
doubling, intermittency and quasiperiodicity routes to chaos [1,2,3,4,5,6]. The 
outcome of this extension, obtained basically from an a priori approach, are 
significant new results that enlarge considerably the previous knowledge. For 
example, it is now clear that all the properties found for the sensitivity to ini- 
tial conditions are associated to the discontinuities of Feigenbaum's universal 
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a scaling function for the bifurcation cascade onset of chaos [5] (with an anal- 
ogous extension for the quasiperiodic critical attractor [6]). In the language of 
Tsallis' functional forms these properties can be read as Tsallis' generalization 
of the Lyapunov coefficient that exhibits a new spectrum of behaviors ranging 
from weaker to stronger than exponential laws [2,3,4]. Also, precise contact 
with the formahsm of Mori's dynamical transitions at the transition point and 
the fixed g-parameter in Tsalhs' formahsm has been estabhshed [5,6]. In view 
of this relation on can conclude that the dynamics for the period doubling 
and the quasiperiodic critical attractors is in fact characterized by an infinite 
number of Mori's dynamical transitions [5,6], in place of the single one known 
previously. In relation to a recent debate in the literature [27], let us then 
conclude that the novel exact knowledge on the critical dynamics obtained in 
[1,2,3,4,5,6] puts into a different fresh context the fundamental scaling laws 
discovered in the 80's for these nonlinear dynamical systems. 

As an example of how to implement the detailed knowledge of the dynamics 
at the period doubling onset of chaos, we have reported [7,8] that specific par- 
allels with the behavior of a glass former close to the glass transition can be 
established if one adds a noise term to the logistic map equation. Of course 
through this very simple one-dimensional system we can aim only at qualita- 
tive analogies with respect to the behavior of a supercooled liquid close to the 
glass transition, nevertheless these similarities have the appeal of allowing a 
precise numerical and analytical inspection. 

This paper is dedicated to Alberto Robledo, with whom I had the great privi- 
lege and pleasure of sharing personal and professional views during long term 
collaborations and engaging discussions. 
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